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Abstract 

This paper provides novel analytic expressions for the incomplete Toronto function, Tg (m, n, r), and 
the incomplete Lipschitz-Hankel Integrals of the modified Bessel function of the first kind, Iem,n{a, z). 
These expressions are expressed in closed-form and are valid for the case that n is an odd multiple of 
1/2, i.e. n±0.5 S N. Capitalizing on these, tight upper and lower bounds are subsequently proposed for 
both Tb(to, n,r) function and Iem,,nia-, z) integrals. Importantly, all new representations are expressed 
in closed-form whilst the proposed bounds are shown to be rather tight. To this effect, they can be 
effectively exploited in various analytical studies related to wireless communication theory. Indicative 
applications include, among others, the performance evaluation of digital communications over fading 
channels and the information-theoretic analysis of multiple-input multiple-output systems. 

Index Terms 

Incomplete Toronto function. Incomplete Lipschitz-Hankel Integrals, Marcum Q-function, upper and 
lower bounds, fading channels. 
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I. Introduction 

Special functions are undoubtedly an inevitable mathematical tool in almost all areas of 
natural sciences and engineering. In the wide field of telecommunications, their use in numerous 
analytical studies often renders possible the derivation of explicit expressions for important per¬ 
formance metrics such as channel capacity and error probability. Furthermore, their corresponding 
computation is typically not considered laborious since the majority of them are included as built- 
in functions in widely known software packages such as Maple, Mathematica and Matlab. 
Based on this, both the algebraic representation and computational realization of any associated 
analytical expressions have been undoubtedly simplified. 

The incomplete Toronto function and the incomplete Lipschitz-Hankel integrals (ILHIs) ap¬ 
pear, among others, in various analytic solutions of problems related to wireless communications. 
They were both proposed a few decades ago and they are denoted as Tsim, n, r) and Zem,n{a-, z), 
respectively 0]. The incomplete Toronto function constitutes a special case of the complete 
Toronto function, which was initially proposed by Hatley in [l2l. It also includes as a special 
case the Marcum Q-function and has been used in studies related to statistics, signal detection and 
estimation, radar systems and error probability analysis, [Q-llSl. Its definition is typically given 
in integral form which involves an arbitrary power term, an exponential term and a modified 
Bessel function of the first kind while alternative representations include two infinite series, 
[|6l . In the same context, the ILHIs belong to a class of incomplete cylindrical functions that 
have been largely encountered in analytical solutions of numerous problems in electromagnetics, 
0, m and the references therein. In communication theory, they have been used in recent 
investigations related with the error rate analysis of MIMO systems under imperfect channel 
state information (CSI) employing adaptive modulation, transmit beamforming and maximal 
ratio combining (MRC), [|9l . 

However, in spite of the evident importance of the TB{rn,n,r) functions and the Zem,nio.,z) 
integrals, they are both neither tabulated, nor included as built-in functions in the aforementioned 
popular software packages. As a consequence, they appear inconvenient to handle both analyt¬ 
ically and computationally. Motivated by this, the aim of this work is the derivation of novel 
analytic results for TB{m,n,r) and ^)- In more details, explicit expressions and upper 
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and lower bounds to the Tsirn, n, r) function and z) are derived for the case 

of n + 0.5 G N. The offered results are expressed in closed-form and have a tractable algebraic 



are derived for the case 


representation which ultimately renders them useful for utilization in various analytical studies 
associated to wireless communications. Indicatively, such studies include, among others, the 
derivation of explicit expressions for important performance measures, such as channel capacity 
and probability of error, in the wide field of digital communications over fading channels and 
the information-theoretic analysis of MIMO systems, among others lfT0l - llT4l . 

The remainder of this paper is structured as follows: Section II revisits the definition and 
basic principles of the Tsim, n, r) function and the Zem,nia, z) integrals. Subsequently, Sections 
III and IV are devoted to the derivation of novel expressions and upper and lower bounds, 
respectively. Finally, discussion on the potential applicability of the offered relationships in 
wireless communications along with closing remarks, are provided in Section V. 

II. Definitions and existing representations 
A. The Incomplete Toronto Function 

The incomplete Toronto function is defined as. 



( 1 ) 


where /„(.) denotes the modifies Bessel function of the first kind and order n. For the special 


case that n = (m — l)/2, it is equivalently expressed in terms of the Marcum Q-function, 
Qm{a,b), as follows 



( 2 ) 


Two alternative representations for the TB{a,b,r) function, in the form of infinite series, were 
reported in O, namely. 



( 3 ) 


and 


'Only the /„(*) function is considered in the present work. 
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where, 


Tsirn^ n,r) = r 


_ ^2(n-a+l)g-r' 


E 

k=0 


r^^'y{a + k, B) 
k\{n + k)\ 


(4) 


n = E 


(n + l)ii\ 


and a = (m + l)/2. 

The notations (a)^ and 7 ( 0 , 0 ;) denote the Poehhammer symbol and the lower ineomplete 
gamma funetion, respeetively, IIT5l - lfT^ . Notably, equations Q and (0) are exaet while, their 
algebraie representation appears to be relatively traetable. Nevertheless, their infinite form raises 
eonvergenee and truneation issues. 


B. The Incomplete Lipschitz-Hankel Integrals 
The general ILHI is defined as, 

Zem,n{(^,z)^ f ( 5 ) 

Jo 

where m, n, a, z may be also eomplex |171, |[8l. The notation Zn{x) denotes one of the eylindrieal 
funetions Jn{x), In{x), Yn{x), Kn{x), Hl{x) or H‘l{x), [[U. An alternative representation for the 
ILHIs of the first-kind modified Bessel functions, was recently reported in dU. This representation 
is given in terms of the Marcum Q-function and is expressed as 

m n+l 

z) = AlJa) + e-” 

i=0 j=0 

+ ^m,n(a)Qi v^\/ a + s/a'^ -1^ (6) 

where the set of coefficients and can be obtained recursively. As already 

mentioned, the above relationship has been shown to be useful in the error rate analysis of 
MIMO systems under imperfect channel state information (CSI). 
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III. An Exact Representation and Bounds for the Incomplete Toronto 


Function 


Recalling Section I, the TB{m,n,r) function is neither expressed in terms of other special 
and/or elementary functions, nor is it included as built-in function in popular mathematical 
software packages. Motivated by this, a novel closed-form expression is derived for the case that 
n is an odd multiple of 1/2. Capitalizing on this expression, novel closed-form upper and lower 
bounds are subsequently deduced. 

A. A Closed-Form Solution for the TB{m,n,r) Function. 

Theorem 1. For m,r,B G n±iGN, m±|GM and m > n, the following relationship 
holds, 




(7) 


where L = m — n — ^ and 7(0,0;) denotes the lower incomplete gamma function IfTI . 


Proof. By setting in x = 2rt and assuming n + ^ G M, the corresponding In{x) function 
can be re-written according to |[l5l eq. (8.467)], namely. 



( 8 ) 


By substituting in ([T]) and making use of the basic identity: (a ±6)^ = ±2a6-|-6^, one obtains 




( 9 ) 
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where L = m — n — k — Evidently, a elosed-form solution to the above expression is subject 
to evaluation of the two involved integrals. To this end, with the aid of [HH eq. (1.3.3.18)], for 
the case that L = 2k with fc G N, one obtains 


pB pB ^ 

/ = V 77 


Ua'- 


-I 


f*B-\-a 


f*B—a 




To this effect, equation (9) can be re-written as follows 


T)^ / x'e + I x'e 

Jo Jo 

( 10 ) 


^ (n + k-l)\L\2 '■ 

TB{m,n,r) = ^ 


(- 1 ) 


m—k—l 


^ot"o {n-k-\)\l\{L-l)\ 

/ B-\-r pB—r 

Je-^^dt - (- 1 )'= / t^e-^^dt 


( 11 ) 


Evidently, the above integrals can be solved in terms of the lower incomplete gamma function 
according to ifTSl eq. (3.381.3)]. Therefore, equation (jV]) is finally deduced and the proof is 
completed. ■ 


B. Upper and Lower Bounds for the Tsirn, n, r) Function 

Novel bounds to the incomplete Toronto function may be straightforwardly deduced from 
Theorem 1. 


Corollary 1: For m,r,B,b G M’*' and m>n, the following inequality holds, 

TB{m,n,r) > Tg (m, [n + 0.5] — 0.5, r) (12) 

where TB{m, [n + 0.5] — 0.5,r) is given in closed-form in ([7]) since it meets the condition 
n + I G M. 

Proof: The incomplete Toronto function is strictly decreasing with respect to n. To this effect, 
for an arbitrary real positive value a, it follows that TB{m, n + a,r) < TB(m, n, r). As a result, 
by recalling that the incomplete Toronto function can be expressed in closed-form for n -f ^ G M, 
by upper ceiling the TB{m,n,r) according to the identity \a + 0.5] — 0.5 > a, lower bounds 


dx 
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the function as in (11) and thus, the proof is completed. ■ 

Corollary 2: For m, r,B,b E M"*" and m >n, the following inequality holds, 

TB{m,n,r) < Tg (m, [n — 0.5J + 0.5, r) (13) 

where Tb (m, [n — 0.5J + 0.5, r) is given in closed-form in ([7]) is given in closed-form in ([7]) 
since it meets the condition n + | G N. 

Proof: The proof follows immediately from Theorem 1 and Corollary 1. ■ 



Fig. 1. Behaviour of the exact solution and the performance bounds to TB{m, n, r) for m = 1.0 and different n 


C. Numerical Results 

The validity of the derived closed-form expression and the general behaviour of the offered 
bounds are illustrated in Figures 1 and 2. Specifically, the behaviour of ([7]) with respect to r, is 
depicted for Figure 1 for n = 0.5 and m = 1.0. Results obtained from numerical integrations 
for m = 1.0 and n = 0.4, n = 0.5 and n = 0.6 are also demonstrated for comparative purposes. 
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In the same eontext, equation (|7]) is depicted in Figure 2 for m = 3.0 and n = 2.5 along with 
numerical results for the cases that n = 2.4, n = 2.5 and 2.6. Evidently, one can observe that 
(|7]) is in exact agreement with the corresponding numerical results while the overall tightness of 
the derived bounds is shown to be quite adequate over the whole range of values of r. 

IV. An Exact Representation and Bounds eor the the Incomplete 

Eipschitz-Hankel Integrals 

Eikewise the TB{m,n,r) function, the lEHIs are neither tabulated, nor are they built-in 
in widely known mathematical software packages. However, their algebraic form constitutes 
possible the derivation of a closed-form expression for the case that n is an odd multiple of 1/2. 


A. A Closed-Form Solution for the Iem,n{z;a) Integrals 

Theorem 2: For m,r,B G M’*', n -f 0.5 G M and m> n, the following closed-form relationship 
holds, 




E 

fe =0 


{n + k-\)\ 


y/HM (n — k — i)! 2^+2 


(a-1)^) , , 1^n+i7(^, + 

(a-l)f' (a+l)'> / 

(14) 


where 


P = m — k + 


1 

2 ‘ 


(15) 


Proof: By expressing the In{,x) function in with its closed-form representation according to 
[fTSl eq. (8.467)] and substituting in Q, it immediately follows that 




-k-i 


[n + k -}^\2 
b; {n-k-\)\ 


k=0 


-1)^ r 

Jo 


x^e-^^^e^dx -f (-1)^+2 


xP^-aXe-xdx 


( 16 ) 
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Fig. 2. Behaviour of the exact solution and the performance bounds to TB{m, n, r) for m = 3.0 and different n 

The involved integrals in (fT^ have the form of the lower ineomplete gamma function. Hence, 
by carrying out some necessary algebraic manipulations and with the aid of [IT5l eq. (3.381.3)], 
one obtains which completes the proof. ■ 

Remark: The present analysis was limited in the consideration of the In{x) function. Nev¬ 
ertheless, similar expressions can be analogously derived for the case of the Bessel functions 
Jn{x), Ynix), Kn{x) as wcll as the Hankel functions, Hn\x) and hL^\x), [[I]|. 

B. Upper and Lower bounds for the z) Integrals 

Corollary 3: For m, r,B,n E M'*' and m >n, the following inequality holds 

z) > 1 (a, 2 :) (17) 

where /^ „_i(a,; 2 ) can be expressed in closed-form according to (13) since it always meet the 
condition n + | G N. 
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Proof: it is noted that the Im,n{.o.,z) integrals are monotonically decreasing with respect to 
n. Thus, for an arbitrary real positive value a, a e M"*", it follows that /m+a,n(a, z) < ^)- 

Thus, by recalling the property fn + 0.5] —0.5 > n the closed-form lower bound in (fTTl) is 
deduced. ■ 


Corollary 4: For m, r,B,n E M'*' and m >n, the following inequality holds 

z) < (a, ( 18 ) 

where + i (a, z) can be expressed in closed-form according to (13) since it always meet 

the condition n -f | G M. 

Proof: The proof follows immediately from Theorem 2 and Corollary 3. ■ 



Fig. 3. Behaviour of the exact solution and the performance bounds to Iem,n{a, z) for different values of n and m 
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C. Numerical Results 

The validity and behaviour of the offered results are demonstrated in Figure 3. One can observe 
the exactness of (fT4l) along with the evident tightness of the proposed bounds. Importantly, the 
achieved tightness holds over the whole range of parametric values. 

V. Conclusion 

In this work, explicit representations and upper and lower bounds were derived for the 
incomplete Toronto function and the incomplete Lipschitz-Hankel integrals of the modified Bessel 
function of the first kind. The offered results are novel and are expressed in closed-form. This is 
sufficiently advantageous since it renders them suitable for application in various studies relating 
to the performance analysis of digital communications over fading channels, among others. 
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